AD-A075  745  ARMY  MISSILE  COMMAND  REDSTONE  ARSENAL  AL  TECHNOLOGY  LAB  F/G  12/1 

DETERMINATION  OF  DIGITAL  CONTROL  SYSTEM  RfSPONSE  BY  CROSS-MULTI— ETC (U) 
may  79  S  M  SELTZER 

UNCLASSIFIED  DRSMI-T-79-58  NL 

of  ^^1^3  HH 


cznvrw 


U.S.  ARfflY 

missne 

RESEARCH 


ROD 


DEVELOPmEAT 

commAAD 


Rcdsiona  Arwnai.  Alabama  36109 


TECHNICAL  M90RT  T-Tt-M 


DETERMINATION  OF  DIGITAL  CONTROL 
SYSTEM  RESPONSE  RY 
CROSS-MULTIPLICATION 


SM  Saltttr 

Guidance  and  Control  Directorate 


2f  May  107* 


0 


D  D  C 

f?f7Di?mnn 


OCT  .10  UTS 


UAsM  U  U 


Appro***  If»  Public  Willi  OsrntHitiOn  Unlimrtad 


79  10  29  143 


» 


DISPOSITION  MSTDUCnOMS 

DESTROY  this  REPORT  WHEN  IT  IS  NO  LONGER  NEEDED  DO  NOT 
RETURN  IT  TO  THC  ORIGINATOR 


ThC  F INOINOS  IN  this  REPORT  ARC  NOT  TO  M  CONSTRUED  AS  AN 
OFFICIAL  DEPARTMENT  OF  THC  ARMY  POSITION  UNLESS  SO 
DCSIGNATEO  BY  OTHCR  AUTHORl/EO  DOCUMENTS 


TRADE  NAMES 

USE  Of  TRADE  NAMES  OR  MANUFACTURERS  IN  THIS  REPORT  DOES 
NOT  CONSTITUTE  AN  OFFICIAL  ENDORSEMENT  OR  APPROVAL  Of  THE 
USE  OF  SUCH  COMMERCIAL  HARDWARE  OR  SOFTWARE 


REPORT  DOCUMENTATION  PACE 


Determination  of  Digit*!  Control  System 
Response  By  Cross-Mul t tplicat lon# 


Approved  (or  public  release 


This  report  describes  the  Cross-Multiplication  technique  (or  obtaining  the 
response  o(  a  digital  control  syatern  fro*  ita  c lowed- loop  transfer  (unction. 
When  the  latter  la  expreaaed  in  the  complex  z -domain,  the  reaponae  ia  obtain¬ 
ed  at  the  sampling  inatanta.  If  it  ta  desired  to  know  the  reaponae  between 
the  sampling  instants,  the  submultiple  method  siav  be  adapted  to  the  cross- 
multiplication  technique.  When  the  transfer  function  ia  expreaaed  in  the 
modified  z-transform  domain,  the  reaponae  la  obtained  at  any  point  within  each 


sampling  instant.  The  three  approac 
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I.  INTRODUCTION 


I  his  report  describes  a  technique  lor  obtaining  l he  response  ol  a  digital  control  system  It  1% 
assumed  that  the  closed-loop  transfer  (unction  is  asailable  in  the  /-  or  modified  /-transform 
domain  The  numerator  and  denominator  ol  each  side  of  the  transfer  function  are 
cross-multiplied  The  Real  Translation  theorem  is  then  applied  to  the  result,  yielding  a 
difference  equation  in  the  time-domain  1 1  his  may  be  solved  for  the  system  response  in  terms 
of  the  reference  (or  other)  inputts)  to  the  system  as  well  as  in  terms  of  system  state  initial 
conditions 


two  different  modifications  to  the  basic  technique  are  described  one  using  the  suhmultiple 
method  and  one  using  the  modified  /-transform  technique  I  hese  are  applied  when  it  is  desired 
to  determine  intra-sampling  responses  of  the  system  All  three  techniques  are  applied  to  a 
single  csamplc  A  summary  of  the  techniques  and  their  application  is  provided  at  the 
conclusion  of  the  report 


2.  RESPONSE  AT  SAMPLING  INSTANTS 


It  is  assumed  that  a  given  digital  or  samplcd-data  system  can  be  described  by  a  closed-loop 
transfer  function  that  relates  the  controlled  output  of  the  system  to  the  reference  input  If  there 
is  more  than  one  input,  (he  technique  can  also  be  applied  to  the  resulting  sum  of  closed-loop 
transfer  functions  relating  the  controlled  output  to  each  of  the  inputs  Although  the  report 
refers  only  to  a  single  controlled  output,  the  technique  can  be  applied  to  find  any  sy  stem  state  if 
it  is  related  to  the  inputs  to  the  system  in  the  /-domain  !  hese  relationships  may  be  derived  by 
using  any  of  the  standard  techniques  (such  as  signal  flow  graphs)  or  by  the  newly  developed 
SAM  (Systematic  Analytical  Mcthodl  technique 

It  is  assumed  that  the  state  whose  response  is  desired  is  denoted  in  the  /-domain  as  (.  (/).  where 

C(z)  y  jc  (t)  |  .  (1) 


I  RC  Kuo.  ina1\u*  anj  S\mhru%  #•/  \mmpir4 l>otn  (  fntrul  .Vn/rwit,  Pfcnlut  Hill  Sm  Jfoo.  I96J 

J  NM  Nclt/er  N  4  U  4n  Ahrtnainr  Samp  Ini-  I'm*  Sifnml  (irapht.  \  S  Arm}  Motile  Research  and  I  >rv  clop  men  1 

(  ommiod  Red*!  on*  Artenal.  Alabama.  I  ethnical  Report  !  Mat  I**7* 


I  he  script  /  denotes  the  operation  ot  taking  the  /-transform  further  utume  that  there  is  onl> 

one  input  into  the  ttticm  Ki/l. 

where 


R  ( z) 


(2) 


The  relationship  between  ('(/>  and  Kt/i  usually  can  be  expressed  as  a  closed-loop  transfer 
function  (or  several  such  transfer  functions!  which  is  a  ratio  of  two  polynomials  in  /.  i  e  . 


C  ( z) 
rTzT 


M 

r 

Jjl9_ 


aj  z‘ 


N 

z  bR «' 

k-0 


(3) 


where  coefficients  a  and  b  represent  the  system  parameters  The  procedure  for  finding  the 
response  at  sampling  instants  by  the  cross-multiplication  method  consists  of  three  steps 
Step  I  C  ross-multipls  the  numerators  with  the  denominators  of  Equation  (3).  yielding  the 
expression. 


£  bk  2  c,2>  * 

k  «0 


y  a^  z^R(z) 

j-o 


(4a) 


b0C(z)  +  b^  zC ( z)  +  ...  +  bkzkC(z)  +  ...  +  bNzNC(z) 


i  M 

aQR(z)  +  a^zRfz)  ♦  ...  +  a^zJR(z)  ♦  ...  +  aMz  R(z)  . 


(4b) 
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Step  2  Each  side  of  Equation  (4a»  or  (4hl  in  divided  b>  bs/s  I  he  mulling  equation  it  then 
solved  lor  the  (.'(/>  term  that  is  nol  multiplied  hv  a  non-zero  power  of  /.  i.e  . 


coi  -  ^  *'Nr<*>  *  ^  *  . 

N  N 


z^~NR(z) 


.  *"  zM-NR(.)  b°  ,-"cl,)  bl  «1_HC(«) 


k  zk~NC(z) 


.  "-1  z-’c.z,  . 


Step  f  Applv  the  Real  Translation  I  hcorcm  to  Equation  (5l.  recalling  that 

DO 

y'1  |c ( z) |  -  cMt)  3-  ^Pc(nT)  6  (t-nT) 


c(t-kT) 


}  ’ 


Z  C  (  Z)  . 


I  he  asterisk  is  used  to  indicate  a  variable  that  has  been  sampled,  and  btt-n  I )  denotes  a  Dirac 
function  occur  ing  at  the  instant  l=n  1  I  he  resulting  value  of  ctnl )  at  each  sampling  instant, 
nl.  then  becomes 

a  a. 

c(nT)  *  c~  r  (  (n-N)  T]  +  ,-±  r((n-N+l)T]  +  ... 

N  dN 

+  r-i  r((n-N+j)TI  •»■  .  .  .  ♦  r-^  r  f  (n-N+M)  T) 


b  b. 

c~  C  [  (n-N)  T]  -  r-i  c  (  (n-N+1 )  T]  -  ... 


D  . 

c  (  (n-N+k) T)  -  ...  +  cI(n-l)T](  <7) 


5 


where  j.  k.  M.  N.  and  n  are  integer* 

Ihe  advantage*  of  the  form  of  Equation  (7)  are  three-fold 

•  1  he  value  of  ctn  I ).  for  any  l  =  nT.mav  he  obtained  for  an>  form  of  rttl.  whether  or  not 
it  is  “/-transformable 

•  Fhe  expression  for  ctn  1  tdoes  not  have  lobe  recalculated  every  time  the  form  of  rtl) 
change*,  a*  i*  the  ca*e  w  hen  the  response  i*  determined  by  the  partial  fraction,  power  series,  or 
inversion  formula  methods'. 

•  I  he  form  of  the  expression  for  ctnf  i  permit*  the  inclusion  of  initial  condition*,  such 
a*  ctO).  if  they  exist 


FICTITIOUS  SAMPLER 
c(r) 


Fiqure  1.  Closed-loop  sampled-data  system. 


Example  I  Coven  the  closed-loop  sampled-data  system  of  f  igure  /.  'the  closed-loop  transfer 
function  easily  is  found  to  be 


C (z)  _  G ( 2 ) 

R  (  z)  “  1  +  HG  (z)  * 


(8) 


If  < it s >  and  H(s)  arc  given  to  be 

G(s)  *  1/  s(s  +  1) 


and 

H  (s)  *  1, 


(9) 
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respectively.  their  /-transforms  arc 


G(z) 


-T 

z  (1  -  e) 

(z  -  1)  (z  -  e"T) 


(10) 


and 


HG ( z)  -  G(Z)  . 


(11) 


II  one  substitutes  the  expressions  of  I  quationsi  I0l  andt  1 1 1  into  I  quationtKi.  one  obtains  the 
closed-loop  (rangier  function. 


C ( z)  _  (1  -  e~T)  z _ 

R  ( z )  z2  -  2e"Tz  ♦  e 


(12) 


One  ma>  nos*  apply  the  Ihree  steps  prescribed  lor  Ihc  cross-multiplication  procedure, 
obtaining 

Step  I 


z  C(z)  - 


2e  TzC(z) 


e“TC(z) 


(l-e_T) 


ZR(Z) 


(13) 


Step  2 


C(z)  -  (1  -  e_T)z'1R(z)  ♦  2e-Tz-1C(z)  -  e‘Tz“2C(z)  •  (14) 


Step  .1 


c(nT)-(l  -  e"T)r  [(n-l)Tl  ♦  2e"Tc  l(n-l)T]  -  e"Tc  l(n-2)T]. 

(15) 

If.  as  in  pp  145-147  of  Kuo’s  book,  rit)  is  assumed  to  be  a  unit  step  input,  the  system  is  assumed 
to  be  initially  at  rest,  and  the  sampling  period  T  is  assumed  to  be  I  second,  application  of 


Equation  (13)  readily  yields  the  following  value*  for  c(nT) 
efT)  ■  0  6321.  the  value  of  c*tt)  at  I  see 
ct 2 T)  *  I  0972  the  value  of  c*(t)  at  2  sec 
C(3T)  *  I  2067.  the  value  of  c*(t)  at  3  sec 

ctnl )  *  0  6321  +  0  736  4(n-l)T|  •  0  *68jc  (n-2)T).  the  value  of  c*(t)  at  t  n  sec 
Theve  value*  correspond  to  those  obtained  h>  more  tedious  means  in  page*  143-147  of  Kuo's 
book. 

3.  RESPONSE  BETWEEN  SAMPLING  INSTANTS  USING  THE 
SUBMULTIPLE  METHOD 


If  it  i»  desired  to  find  the  intra-sampling  response  of  the  same  type  digital  system  described 
in  Section  2.  it  ma>  be  accomplished  bv  applying  the  submultiple  method  found  in  page*  X3-K6 
of  Kuo’s  book  Briefly,  let  c  *T  represent  the  value  of  the  response  ct 1 1  at  the  instant,  t  * 
^7  .  where  n-l  represents  the  number  ol  intrasampling  response*  desired  in  nan  integer 

with  value  greater  than  unity)  If  m  is  an  integer,  the  sampling  period  within  which  the 
submultiplcs  arc  to  be  determined  isdenoted  a*  ml  I  he /-transform  of  c  yj  maybefound 
from  the  ordinary /-transform  in  the  following  manner  In  essence,  c  is  the  output  of  a 

fictitious  sampler  which  sample*  n  time*  as  fast  as  the  real  sampler  The  /-transform  of  that 
output  is  defined  as 


y  {“s1’}*'  c(z>n 


C(z) 


z-z 

r 

T-Tr 


w  here 


Zn  *  Z 


1/n 


and 


T  «  T/n. 
n 


(16) 


(17a) 


(17b) 


Now  the  closed-loop  expression  of  I  quation  (3)  may  be  altered  to  read 


C(z,n 
R  (z) 


k 


N 

X 

k-0 


bk  zn 


(18) 
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I  he  submulnple  modification  to  the  bane  method  also  consists  of  three  steps 


Step  I  Cross-mult  iply  the  numerator  and  denominators  of  tqualion  (181.  yielding  the 
expression. 


N 


k-0 


VnC(*> 


n  "  E  V* 
j-0 


R(z) 


(19a) 


or 


b0C(z,n  *  bl'nCII'n  4 


N, 


Vn  C|I,„4"-  4  Vn  C(I,n 


-  a0B(z)  +  alznR(2)  + 


4  M 

♦  a  z  JR(z)  ♦  ...  ♦  aMz  R(z)  . 
j  n  pi  n 


(19b) 

Step  2  I  aeh  nde  of  I  quahon  ( I  S»a  >  or<  I9b>  is  divided  by  bs/'  I  he  resulting  equation  is  then 
solved  for  (he  C(/>.  term  that  is  not  multiplied  by  a  non-zero  po*er  of  /„.  i  c  . 


C  ( Z ) 


-  —  z‘N  R  ( Z)  +  z  •(N’l)R(Z)  ♦  ... 
bN  n 


n  bo  n 


4  il*  -(N"j,R(Z>  ♦  ... 
bN  n 


M  ,  - (N-M) 

3n 


+  p  zn  * . R(Z)  -  r-5.  z  N  C(z) 

n  bN  n  n 


ri  z  ’(N-1,c(z,n 

bN  n 


k  Z  “(N“k)C(z) 
N  " 


z  _1C ( z)  .  (20) 

N  n 
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Step  3  Similar  10  the  procedure  ol  Section  2.  the  value  of  an  at  the  n“  tubmultiple  ol  ihc 
sampling  instant  ml  is 

[<=^1  ♦  ••• 

. 

where  all  values  ol  rf  |  are  equal  to  zero  except  lor  those  values  of  r  at  integral  multiples  of  I. 


OV  I  met  .nr 


r(»T) 


t-*T  for  |  -  0.  I.  2. 


i  integers! 


ami  j.  k.  M.  m.  N.  n  are  integers 

E  sample  2  The  same  sampled  data  system  is  used  in  this  example  as  in  Example  I  It  may  he 
shown,  for  the  example  at  hand,  that 

C(*)n  •  G(Z) n  E ( z) ,  (23) 

where  (»</».  is  found  in  the  same  manner  as  demonstrated  in  Equation!  I6i  Erom  hgurr  /  it  is 
seen  that 


E(Z)  - 


R(Z) 

4  HG ( Z 1 


Substitution  of  Equation  <24l  into  Equation  (23l  leads  to 


C(z)  G ( Z) 

_ n  _  n 


1  4  HG ( Z ) 


Replacing  <  by  <•  and  T  by  T.  in  Equation  (10)  trade  to  the  eapreuion. 


V1  •  •~T|'» 

(,n  '  11  <‘n  *  e'Tn' 


If  one  substitutes  Equations  4 1 7)  into  filiation  (III  and  substitutes  the  mulling  equation  and 
Equation  (26)  into  Equation  (23).  one  obtain*  the  c lined -loop  cupicuion. 

C  ( z)  z  ( 1  -  e-Tn)  (x  -  1)  (x  -  e'T) 

n  n  _  _  <111 


R(zi  (x  -1) (z  -e”^n) (z 

n  n  n 


x  2n-2e"Tz_n*e"T) 


I  he  denominatot  o(  this  expression  is  of  course  (when  equated  to  /eroi  the  characteristic 
equation  It  mas  he  expanded  into  the  following  polynomial  in  u 

s  2n*2  . ,  ,  -T  ,  2n+l  .  -T  2n 

D(z  )  -  x„  -  (1  ♦  e  n)  z  ♦  e  n  z 

n  n  n  n 


-  „n‘»  *  2.-t  (l  ♦  "♦*.  2.-,T*Tn'znn 


♦  <="T»n2  -  *'T(1  ♦  *'*")»„  ♦  e'(T*T"’. 


Similarly,  the  numerator  of  I  quation  (27|  mas  he  expressed  as  the  polynomial 


N(zn)  -  (1  -  e*Tn)  [zn2n+1  -  (1  ♦  e"T)znn+1  ♦  e"Tzi 


Step  I  If  Equation*  |2K)  and  (24)  arr  *uh*tituted  into  Fquation  (27)  and  the  mulling 
numerator*  and  denominator*  crou-tnultiplied.  one  obtain*  the  equivalent  of  Equation  ( 19b): 

z  2n*2C(z)  -  (1  ♦  e"Tn)  z  2n+1C(z)  ♦  e-Tnz  2nC( z) 

n  n  n  n  n  n 

-  2e"Tz  n*2C(z)  ♦  2e"T(l  ♦  e‘Tn) zn"+1C ( z) 

n  n  n  n 

-2.  -<T*Tn>,  "CL>„ 
n  n 

♦  e”Tz  2C ( z)  -  e‘T(l  ♦  e~Tn) z  C(z)  ♦  e‘ (T+Tn> C ( z) 
n  n  n  n  n 

-  (1  -  e"Tn>  [zn2n+1R(z)  -  (1  ♦  e~T) z  n+1R ( z) 
n  n 


♦  e  Tzn  R ( z) J  . 


(30) 


Step  2  If  one  applie*  the  procedure  of  Step  2  to  Equation  (30)  one  obtain* 

C ( z)  -  (1  -  e"Tn)  It’1  R ( z )  -  (1  ♦  e"T)z  ■(n+1,R(z) 
n  n  n 

<  e_Tz‘(2n*1)  R  ( z )  ]  ♦  (1  ♦  e*Tn)  z  _1C  ( z) 

n  n  n 

-  e"Tnz  “2C(z)  ♦  2e'Tz“nC(z) 

n  n  nr 

-  2e-TU.e-T")In-|n*1>c(.)n.  2.' lT*Tn’*n' 


-  e‘Tz  ” 2nC ( z)  ♦  e_T(l  ♦  e"Tn)2  "(2n+1) 
n  n 

- (T+T  )  _  - ( 2N+2) 


-  -  .  -  C  (z) 

n  n  n  n 


n  z 


C  (z) 


n  . 


(31) 
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Sicp  .V  Finally,  one  may  apply  the  procedure  of  Step  .1  to  obtain  the  time-domain  difference 
equation  that  yields  the  value  of  c*(t)  at  the  instant  of  time,  t»  — -  ,  where  m  is  any  desired 
integer 


■\gain.  as  in  Kuo's  book  and  in  I  sample  I .  it  is  assumed  the  system  is  at  rest  initially,  and  the 
sampling  period  I  is  one  second  One  may  desire  to  know  two  intra-sampling  values  of  the 
output  c*tti  In  that  case,  n  2+1  *  If  one  substitutes  these  numerical  values  into  Equation 

(321.  one  obtains  the  equation. 

c(j)«0.2835r(^)-0.4866  r(^) 

+  0.1043r(^)+  1.7165  c(^) 

-0. 7l65c(^)+0. 5272c(?^)-l.  2630c(^) 
+0.5272c(^) 

-0.3679c(^)+0.6315c(^)  -0.2636c(?^)  . 


(33) 


Il  is  worthwhile  to  pause  a  moment  and  consider  the  power  of  the  difference  equation  just 
obtained  First,  the  output  c  [^’  can  be  obtained  for  any  deterministic  reference  input  rft). 
It  only  need  be  specified  at  the  sampling  instants  Second,  initial  conditions  or  instantaneous 
changes  can  be  accommodated  readily  by  the  equation  Third,  difference  equations  are 
particularly  amenable  to  programming  on  desktop  calculators  or  digital  computers  I  hese 
three  adsantages  are  not  enjoyed  by  all  methods  found  in  the  standard  textbooks  for 
determining  the  statets)  of  a  sampled-data  system.  The  submultiplc  method  docs  suffer  from 
the  drawback  of  hasing  to  specify  beforehand  the  number  n-l  of  the  inlrasampling  instants  for 
which  it  is  desired  to  know  a  given  state,  such  as  eft)  in  this  case 

In  the  case  at  hand,  a  brief  demonstration  of  the  calculation  of  eft)  at  instants  t=0.  7  7. 

2T  3.  F .  is  provided  below  A  unit  step  input  is  assumed  so  that  the  results  may  be 

compared  to  those  of  Kuo  in  pages  195-  198  of  his  book  Thesymbolc  •  denotes  the  value 
of  eft)  at  t=  sec 

m=0:  flO)  *  0 

nv=l  0l/3>  -  0.2833  r(0>  +  1.7163  00)  «  0  2835 

m=2:  ct  2  3)  =  0  2835  rf  I  3)  +  I  7165  0  I  3)  =  0  4866 

m=3:  0l)  =*  0  2835  r<2  3)  +  I  7165  02  3)  -  0.7165  0l  3)  *  0  6321 

m=4  04  3)  =  0  2835  rt  1 1  -  0  4866  r<0)  +  I  7165  0  I  >- 0  7165  0  2/3)  +  0  5272  0  I  /  3)  =  0  8407 

m=5  05  3)  *  0  2835  ri  4  31-04866  rfl  3)  +  I  716504  3)  -  07165  01)  +  0  5272  02/3) 

-I  2630  0  I  3)  =  0  9901 

m=6  0  2)  =  I  0972 
m=7:  07  3)  =  I  1464 
m=8:  08  3)  *  I  1816 
m=9  03)  =  I  2067 


4.  RESPONSE  BETWEEN  SAMPLING  INSTANTS  USING  MODIFIED 
Z-TRANSFORM  METHOD 


If  it  is  desired  to  find  the  mtra-sampling  response  of  the  same  type  of  digital  system 
described  in  Sertion  2.  it  also  may  be  accomplished  by  applying  the  modified  r-transform 
method  In  this  method  it  is  necessary  first  to  determine  the  modified  /-transform  equivalent 
of  the  closed-loop  transfer  function  of  Equation  (3).  The  reader's  mind  will  be  refreshed  (as 
needed) 
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1  ct  eft)  denote  the  response  of  a  digital  system  and  c*(t )  its  sampled  output.  The  value  of  c<t) 
at  the  instant  of  time.  t"(n-A)T,  is  the  value  of  eft)  that  has  been  delayed  by  an  increment  of 
time  AT  The  latter  is  represented  symbolically  as  ctl-AT).  If  the  relation. 


m  *  1  -  A»  0  £  m  £  1  •  (34) 

is  used,  the  value  of  ett)  delayed  by  an  amount  AT  after  the  sampling  instant.  t=nT.  may  be 
denoted  as  c  |l  n- I  +m)T).  The  /-transform  of  c((n- 1  +m)T|  is  termed  the  modified  /-transform 

(  y  m  1  of  **>• ' e  > 

11  {c  (t-T)  )  =  >  ( c  l  (n-l+iti)^)}  « 

(35) 

^m(c  ( t)  }  C(z,m)  -  C  ( z ,  A)  l_JD  . 


It  also  may  be  shown  that 


C(z,m)  =  z  *  ^{c(t+mT)}  ■  z  *  ]  c  ((m+n)Tl  z  n.  (36) 

k«0 


Alternately.  (.(/.mi  may  be  derived  through  its  complex  convolution  definition 


C ( z ,m) 


*^(c(t-AT)x6T(t)  )  |  z=eTs 
-  (yic<t-AT))*«:UT(t)>]  |j.eTs 

-ik  T  c(tic-«t[i-.-*‘-f]-»at 

c-i® 


(37) 


z=e 


Ts 


9 


where  Ai(t)  represents  a  train  of  impulses  each  with  a  unit  area,  the  star  symbol  (*)  denotes  the 
s  complex  convolution  operation.^  denotes  the  operation  of  taking  the  l  a  place  transform,  c 

denotes  the  abscissa  of  convergence  (a  positive  real  number),  and  £  represents  a  dummy 
complex  variable. 

Tquations  (35)  -  (37)  represent  twfo  methods  of  determining  analytically  the  modified  /- 
transform  (it  is  recommended  that  one  look  it  up  in  a  modified  /-transform  table,  if  one  is 
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available).  The  above  relationships  may  also  be  used  to  determine  the  modified  /-transform  of 
the  transfer  function.  G(/.  m) 


For  the  purpose  of  exposition,  assume  the  modified  /-transform  equivalent  of  the  closed- 
loop  transfer  function  is  known  and  is 


M 

£  a.  z^ 
C(z,m)  j„0  3 

“RTzT  "  N 

£  b*  zk 
k»o 


(38) 


where  a,  and  b  may  be  functions  of  m.  The  cross-multiplication  technique  is  similar  to  those 
two  previously  described 


Step  I  Cross-mult iply  the  numerators  and  denominators  of  Equation  (38).  yielding  the 
expression. 


N  M 

23  zk  C(z,m)  -  J3  a^z^R(z) 

3-0 


(39a) 


k-0 


bQC ( z ,m)  +  b1zC(z,m)  +  b2z2C(z,m)  ♦  ...  +  bkzkC(z,m) 
♦  ...  +  bNzNC(z,m)  *  aQR(z)  +  a^zR(z)  +  ...  ♦  a^z^R(z) 


(39b) 


Step  2  Each  side  of  Fquation  (39a)  or  (39b)  is  divided  by  b\/s.  The  resulting  equation  is 
then  solved  for  the  C(/.  m)  term  that  is  not  multiplied  by  a  non-zero  of  /.  i.e., 

C(z.m)  *  ^z"NR(z)  ♦  ^z'(N‘UR(t)+  •••  ♦  ^*"(N‘j)R(z) 


♦  z'(n_m)r(z)  - 


N 


z  (z , m) 
°N 
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% 


v"  Clz.m)  * 


z 


' C ( z ,m)  ♦ 


♦ 


♦ 


b  z_1C(z,m) 
_N-1 

bN 


(40) 


Step  3  the  inverse  /-transform  of  each  term  of  Equation  <40>  is  determined,  using  the 
inverse  of  Equation  ( 39bl. 


c  [(n-l)T.m]  =  c[(n-l*«)Tj-  g-  r[(n-N)T]  +  ^  r[(n-N*l)T]* 

a .  a 

r  £(n-N*  j ) +  .  ..♦  g-  r[(n-N4m)T] 

N  N 


...  ♦ 


■I 


(41) 


b0 

JT  c| 
bN 


£(n-l*m-N) tJ*  gi  c  £(n^»-N)  tJ ♦  ... 
N 

*  £ 


■  l*m-N*k) T 


•**  *  c  [(n-2«J»)T  J  j 


v» here  ).  K.  m.  M.  N  are  integers 


F  sample  t  Again,  the  same  sampled  data  svstem  is  used  in  this  esamplc  as  in  I  samples  I 
and  2  l  he  modified  /  transform  method  .ill  be  applied  A  fictitious  lime  dclas.  e-*'-*' '.  is 
placed  following  the  forward  loop  gain.  (its),  and  a  fictitious  time  adsancc  lof  the  same 
magnitude  as  the  time  dclas )  is  placed  lust  before  the  feedback  loop  gain  (Its)  I  he  cloved  loop 
transfer  function  that  results  from  these  two  addition  elements  being  added  is 


C(z,m)  _  G (z ,m) 

“rTzT  1  ♦  HG(*r 

(  sing  I  quation  i*>i  for  (its),  one  mav  obtain  (it/  mi 


(42) 


-  ?«{rTirVTT> 

.  ( l~e~mT)  z  ♦  (e~tnT-e~T ) 
(z-l)(z-e_t) 
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G  ( z  ,ro) 


(43) 


The  value  for  Hli(f)  is  given  in  Fquation  (9)  •  (10)  which,  with  Equation  (4)).  may  be 
substituted  into  Equation  (42)  to  yield 


C  ( z ,m) 
R(z) 


(l-e*mT)z  ♦ 

(z-1) (z-e"T) 

fi  ♦  _ 

(z-1) (z-e”T) 

.  ♦  (e~aT  -  e~T) 

z2  -  J.-T,  .  .-T 


(44) 


It  may  be  observed  that  Fquation  (44)  is  in  the  vame  form  as  Equation  (J8)  Hence  the 
operations  denoted  as  Steps  I.  2.  and  J  (above)  may  be  applied 

Step  I  The  numerators  and  denominators  of  Fquation  (44)  are  cross-multiplied,  yielding 

z2  C ( z ,m)  -  2e“T  z  C(z,m)  ♦  e“T  C(z,m) 


-  (l-e""T) z  R(i)  ♦  (e‘®T  -  e'T)z"2R(z) . 


(45) 


Step  2  Each  term  of  f  quation  (45)  iv  divided  by  /  .  and  the  resulting  equation  is  rearranged 
to  solve  for  C(/.  m) 


C  ( z .  m )  -  (l-e^ir^U)  ♦  (e-mT-e"T)  z“2R ( z) 


♦  2e_Tz”1C(z,m)-  e~Tz”2 C ( z  ,m)  . 


(46) 


Step  V  Equation  (46)  is  now  transformed  into  the  time  domain,  resulting  in  the  following 
difference  equation 

c  I  (n-l+m)  T  )  *  ( l-e”mT)  rl(n-l)T]  ♦  (e*mT-e'T)  r|(n-2)T] 


♦  2e"T  c((n-2-HB)Tl-e"T  c  ((n-3+m)T  |  •  (47) 
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During  the  n  "'  sampling  period,  i.e  .  (n-l  )T  <l  <  nT,  where  n  is  an  integer,  one  lets  m  assume 
values  between  zero  and  one  to  find  the  inlrasampling  values  of  ctt).  For  instanced  it  is  desired 
to  verify  the  values  obtained  in  Example  2.  let  m=l/3  and  2/3  and  vary  the  integer  n 
incrementally,  starting  at  n=0  Again,  for  simplicity,  let  T  =  I  second 

n  -  0:  c(m-l)  •  0 

n-l:  c(m)  -  (1-e  m)  r(0)  -  1-e  m 

m  -  0  (check  case):  c(0)  -  0 
m  -  1/3:  c (1/3)  -  1  -  0.7165  -  0.2835 
m  -  2/3:  c(2/3)  -  0.4866 
m  »  1:  c(l)  -  0.6321 

n  -  2:  c (1+m)  -  (l-e*ro)  r(l)  +  (e“m-e’1)  r(0) 

+  2e  1  c(m) 

m  -  0  (check  case):  c  < 1 )  -  0.6321 

m  -  1/3:  c { 4/3)  -  (l-e"1/3)  +  (e"1/3-e_1) 

2e-1  c  ( 1/3) 

-  0.8407 

m  -  2/3:  c(5/3)  -  (l-e~2/3)  +  (e*2/3-e_1) 

♦  2e_1  c(2/3) 

-  0.9901 

m  -  1:  c(2)  -  (1-e*1)  ♦  (e_1-e“1)  ♦  2e_1 
-  1.0972 
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Several  anal) deal  techniques  for  obtaining  the  response  of  a  digital  control  system  have 
been  described  Ihey  arc  based  on  a  single  principle  cross-multiplication  followed  by 
applications  of  the  real  translation  theorems  Fach  is  applied  to  a  single  example.  As  a  starting 
point  for  application  of  each  of  the  techniques,  it  is  required  that  the  dynamics  of  the  digital 
control  system  be  described  in  the  /-  or  modified  /-domain. 
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The  advantage*  of  the  three  technique*  over  extant  classical  methods  arc. 

•  The  response  may  he  obtained  for  any  deterministic  reference  input  into  the  system  as 
long  as  its  value  is  known  at  the  sampling  instants  It  need  not  be  described  by  a  differential 
equation,  and  the  /-transform  for  a  specific  reference  input  need  not  be  determined  before 
obtaining  an  expression  for  the  response 

•  Initial  conditions  and  instantaneous  changes  can  be  accommodated  readily  by  the 
equations  obtained  through  use  of  the  cross-multiplication  methods 

•  I  he  difference  equations  obtained  are  particularly  amenable  to  programming  on  a 
desktop  calculator  or  digital  computer 

•  A  detailed  knowledge  of  the  theory  underlying  digital  or  samplcd-data  control  systems 
is  not  required  (although  it  certainly  is  helpful)  by  the  analyst  in  order  to  apply  the  recipes 
described  herein 
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